AC losses in type-II superconductors induced by 
nonuniform fluctuations of external magnetic field 
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Abstract 

Magnetic field fluctuations are inevitable in practical applications of superconductors and it is often necessary 
to estimate the AC losses these fluctuations induce. If the fluctuation wavelength is greater than the size of a 
superconductor, known estimates for an alternating uniform external magnetic field can be employed. Here we 
consider the opposite case and analyze, using a model critical-state problem, penetration of spatially nonuniform 
fluctuations into type-II superconductors. Numerical simulation is based on a variational formulation of the Bean 
model. The analytical solutions, found in a weak penetration limit, are used to evaluate AC losses for two types of 
fluctuations: the running and standing waves. It is shown that for spatially nonuniform fluctuations the losses are 
better characterized by the fluctuation penetration depth than by the fluctuation amplitude. The results can be used 
to estimate the AC losses in flywheels, electric motors, magnetic shields, etc. 

Index Terms 

Hard superconductors. Bean model, AC losses, penetration depth, nonuniform field fluctuations, asymptotic 
solution. 

I. Introduction 



The possible range of currents and magnetic fields as well as the economic gains of implementation of type-II 
superconductors in power transmission lines, current leads, fault current limiters, magnetic shields, bearings, etc. 
Q are often limited by AC losses and the necessity to remove the generated heat out of the system. Thus, application 
of bulk high-Tc superconductors in flywheel systems and magnetic bearings is promising because of no friction 
I between the moving parts and, hence, no energy losses by friction. However, rotating permanent magnets used in 
^ ' such devices always produce somewhat irregular magnetic field; moving field irregularities cause hysteretic losses 
^ and relaxation of levitation property [1]. 

T— I ^ The mathematical models, used to analyze magnetization of type-II superconductors and to evaluate AC losses, 
^ ■ involve highly nonlinear partial differential equations that have been solved mostly for superconductors in a uniform 
Q alternating external magnetic field (see [2]-[6] and the references therein). In many practical situations, however, 
the external magnetic field is not exactly uniform and can be better presented as a superposition of a uniform 
part and spatiotemporal fluctuations with the characteristic length scale less than the size of a superconductor The 
fluctuations are often stochastic but, for example, in the case of a flywheel system, are induced by rotation of 
^ permanent magnets and may be approximated by a running wave. If there are no moving parts as, e.g., in the case 
of magnetic shields or transformers, nonuniform magnetic field fluctuations in the form of a standing wave can, 
probably, serve as a better approximation. 
Q I Our aim is to study the penetration of nonuniform magnetic field fluctuations into a hard superconductor and 
O to evaluate the accompanying AC losses. We start with a convenient for numerical simulations variational general 
I> formulation of the Bean critical-state model (section |nll, then consider the simplest geometric configuration, a 
^ superconductive slab placed between two parallel sheets of external current. Even in this simplest case the problem 
^ ■ becomes nontrivial if the external current and, hence, also the external magnetic field, are nonuniform. We assume 
, the external field fluctuations are induced by a given current in the form of either a running or a standing wave and 
solve, first, the magnetization problems numerically (section HiH) . Simple physical arguments allow us to find also 
asymptotic analytical solutions for small fluctuations (section Hvt . We further extend these solutions by presenting 

L. Prigozhin is with the Department of Solar Energy and Environmental Physics, Blaustein Institute for Desert Research, Ben Gurion 
University of the Negev, Sede Boqer Campus, 48990 Israel (e-mail: leonid@cs.bgu.ac.il) 

V. Sokolovsky is with the Physics Department, Ben Gurion University of the Negev, Beer Sheva, 84105 Israel (e-mail: 
sokolovv @ bgumail.bgu.ac .il) 



o 



X3 



them as the zero-order terms of consistent asymptotic expansions, find the first order corrections (appendices I and 
II), and, finally, determine the asymptotic AC losses for small fluctuations (section 0. We analyze the dependance 
of the leading AC loss term on the first order correction to current density distribution. 



II. Variational formulation of the Bean model 

Let a superconductor occupying domain i7 C be placed into magnetic field He(r, t) induced by a given 
external current with the density Je{r,t) (here r = {x,y,z} and Je = in Q). In accordance with the Faraday 
law, an alternating magnetic flux induces electric field and, hence, an eddy current inside the superconductor. In an 
ordinary conductor, the vectors of the electric field and current density are related by the linear Ohm law. Type-II 
superconductors are, instead, characterized in the Bean critical state model [7] by a highly nonlinear" cunent-voltage 
relation which gives rise to a free boundary problem. 

Let us assume 

E = pj (1) 

and employ the usual Bean model relations determining the effective resistivity of a superconductor, p(r,t), 
implicitly. Namely, the Bean model states that the effective resistivity is nonnegative, 

p{r,t)>0, (2) 

the current density cannot exceed some critical value, 

\J{r,t)\<J„ (3) 

and, if the current density is subcritical, the resistivity is zero: 

|J(r,t)| < J,^p(r,i) =0. (4) 

Since no current is supposed to be fed into the superconductor by an electric contact, the current density inside Q 
should satisfy the zero divergence condition and have a zero normal component at the domain boundary dQ: 

V • J = in f], J„ = on (917. (5) 

To derive a variational formulation of the magnetization problem, we define the set K of possible current densities, 



J GK 




and express the electric field via the vector and scalar magnetic potentials, 

E = -dtA - Vip. (6) 

We further eliminate the scalar potential by multiplying Q by $ — J, integrating over U, and making use of the 
zero divergence condition for functions from the set K, 

(E, * - J) = - {dtA, * - J) , (7) 

for J, $ G iT. Here (u, v) = u • v dil is the scalar product of two functions. 

Since E and J satisfy the current-voltage relations of the Bean model, the scalar product on the left side of (Q 
is nonpositive for any $ G i^. Indeed, using Q-© we obtain: 

E- J = |E||J| = |E|Jc > |E||*| > E • (8) 

Up to the gradient of a scalar function, determined by the gauge and also eliminated by the scalar product with 
$ — J, the vector potential is a convolution of the Green function of Laplace equation, G{r) = l/47r|r|, and the 
total current density J{r,t) + Je(r, t): 

A = /XoG * ( J + Je) 



(it is assumed that the magnetic permeability of superconductor is equal to that of vacuum, /xq-) We arrived at the 
following variational formulation of the magnetization problem: 

Find J G K such that 
(G*at{J + Je},*- J) >0 for any ^eiC, (9) 
J|t=o = JoW, 

where Jq G is a given initial distribution of the current density. 

The problem ^ is an evolutionary variational inequality with a nonlocal operator and has a unique solution 
[8]. This inequality is written for the induced current density alone: the effective resistivity p has been used 
to derive the inequality ^ and then excluded. It has been shown [8] that, in the Bean model, p{r,t) is the 
Lagrange multiplier related to the current density constraint Of course, the same inequality may be written as 
(g * dtJ + -^dtAe, * - j) > for any ^ G K, where Ag is the external vector potential; such formulation can 
be more convenient in some cases. 

We note that, although the variational formulations where the solution is sought as an extremal point of some 
functional are much more familiar, the variational inequalities do appear in many problems of mechanics and 
physics containing a unilateral constraint or a non-smooth constitutive relation (see [9]). The methods for numerical 
solution of variational inequalities are well developed [10]. Nevertheless, solution of Q in the general 3d case is 
certainly a challenging problem not only because of huge number of finite elements needed: the nonlocal operator 
of convolution leads to a dense matrix of the discretized problem, and the zero divergence condition ^ is not so 
easy to account for numerically, although the edge finite elements and tree-cotree decomposition can be helpful (see 
[11], [12]). Similar difficulties are typical also of other formulations of the magnetization problem [2], [13]. That is 
why the magnetization problems were solved mainly for one and two-dimensional configurations. Below, we use a 
model 2d problem to simulate the penetration of nonuniform magnetic field fluctuations into a hard superconductor 
numerically and to calculate the AC losses asymptotically for weak penetration. 



III. Two TYPES OF FIELD FLUCTUATIONS 

To simulate the penetration of magnetic field fluctuations into a superconductor, let us consider the simplest 
possible configuration: a superconductive slab —l<x<l placed between two sheets of external current at 
X = iba, a > I (Fig. l).We assume the external current is parallel to 2;-axis, Jg = Je^z- Having in mind, e.g.. 
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Fig. I. Infinite slab between two parallel sheets of current. 



a superconductor moving along the axis of a solenoid with slightly nonuniform winding, we present the external 
current Jg as a sum of its uniform component /e and a wave moving along the y-axis with the velocity v: 

Jg = [Ig + i^sui{k{y - vt)]] [5{x + a) - 8{x - a)]. (10) 

We will also consider the standing wave of magnetic field induced by the current 

Jg = [/g + ig sin(A:y) sin(27r/t)] [5{x + a) - 5{x - a)]. (11) 



If the sheet current density were uniform, the Bean model equations for slab magnetization could be solved easily. 
It is more difficult to analyze the magnetization in a nonuniform magnetic field induced by the external currents 
(flOt or (II It . and we will use a 2d reformulation of the variational inequality ^ and the numerical method proposed 
in [4] for cylinders of arbitrary cross sections in a perpendiculai^ magnetic field. The method can be employed for 
any distribution of external currents. Taking the cylinder cross section to be a rectangle, —l<x<l, L < y < L, 
we model the development near the slab surfaces of critical current zones that shield the fluctuating component of 
external magnetic field. 

The shielding current inside the superconductor is directed along the z-axis and does not depend on z, J = 
J{x,y,t)ez, so the conditions ^ are satisfied automatically. We may redefine 17 to be the cylinder cross section 
and solve a 2d problem in this domain. One should, however, be cautious: the eddy current running in the positive 
z-axis direction has to return back (no transport current is applied; see [4] for solution of problems with transport 
current). As a trace of the 3d zero-divergence conditions the condition 



/ J{x,y,t)dn = 
Jn 



has to be satisfied for all t. The set of admissible (scalar) current densities becomes 



K=\ <^{x,y,t) 



|^>| < Jc inn, ^dn = o 



n 

The variational inequality Q can also be rewritten for scalar current densities. 

Find J ^ K such that 

(G * dt{J + Je}, - J) > for any ^ e K, (12) 

J\t=o = Mr), 

where G{r) = —-^ In |r| is now the Green function of Laplace equation in and r = {x, y}. 

To solve the problem numerically, we use first the finite difference discretization in time and obtain, for each 
time layer, the stationary variational inequalities 

Find J e K such that {G * {J + Je - J - Je},*^ - J) > for any ^ e K, 

where J, Jg are the values from the previous time layer. It can be shown that these variational inequalities are 
equivalent to the optimization problems 

mill J, J) + (5, J), 

where g = G * {Je — J — Je}- The convolution G * Je can be calculated analytically: 

1 i 

—Ae = G*Je = -leX - ^c"'"' sm{k(y - vt)} smh(kx) (14) 

—Ae = G*Je = -leX - ^-^e'^" sin(27r/t) sin(%) sinh(A;x) (15) 

for currents dlOt and (IIH . respectively, and —a < x < a. 

We finally discretize the optimization problems il3l in space by means of the piecewise constant finite elements, 
take care of the integral constraint J dO, = using the Lagrange multipliers technique, and solve the resulting 
finite-dimensional optimization problems with the remaining constraints by a point relaxation method (see [4], [14] 
for the implementation details). 

Let us assume only an alternating part of the external current is present and le = 0. A permanent field would 
not change the picture qualitatively and, although in applications the stationary component is often much stronger 
than fluctuations', the constant part causes no energy losses. We present the results of two simulations. In both 
cases we assume the superconductor is initially in the virgin state, Jq = 0. 



'Note that strong permanent component makes the relative changes of magnetic field small. This justifies the use of field-independent 
critical current density for simulating the penetration of fluctuations. Otherwise, the Kim model [15] or its modification can be employed. 
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Fig. 2. Running wave. Penetration of the magnetic field induced by an external current Je = is{t) sm{k{y — vt)}[S{x — a) — S{x + a)]ez 
into a cylindrical superconductor with the cross section Q — {—I < x < I, ~L <y< L}. The wave amplitude ie grows linearly from zero 
at t = to its maximal value i^o at t = to and then remains constant. The physical parameters: I = 1/k, L — 10/k, a = 1.1/k, to = 
4:/kv, ieo ~ Jc/k. Shown for t = [0.8, 1.6, 3.2, 5.6, 16.0]/kv are: the boundaries of plus- and minus critical current zones, the magnetic 
field lines, and the cylinder cross section. The parameters of numerical scheme: 50 x 500 regular mesh of piecewise constant finite elements, 
time step 0.8/A;ii. 

In our first example (Fig. ^ the external current is given as a running wave, Je = ie{t) sm{k{y — vt)}[6{x — 
a) — 6{x + a)], with the amplitude ie{t) growing from zero to its maximal value and then remaining constant. The 
magnetic field penetrates from the surface of the superconductor where alternating domains of plus- and minus 
critical current densities appear and start to follow the wave. The shape of these domains stabilizes and, after an 
initial transient period, they completely occupy a near-surface zone of a constant depth and move through this zone 
with the wave velocity v. 

The second example (Fig. |3ll illustrates another typical situation: here the magnetic field fluctuations are induced 
by the external current in the form of a standing wave, Jg = ig sin(/cy) sin(27r/t)[J(x — a) — 6{x + a)]. At 
t = domains of plus- and minus critical current densities, shielding the external magnetic field, appear at the 
superconductor surface and start to propagate inside. When the external field reaches its maximal strength att = r/4, 
the propagation stops (here T = 1//). As the field becomes weaker, the boundary of the critical current regions 
does not, however, retreat. Instead, similarly to the case of an alternating uniform external field, to compensate 
the decreasing external field there appear surface domains of the opposite critical current densities. These domains 
propagate inside, sweep out the previous ones at t = ST/A, and the process becomes periodic. 

IV. Asymptotic solution 

Let an infinite slab — / < x < / be placed into a magnetic field produced by the sheet current (fTOt or (fTTl . We 
want to find the established periodic in time distribution of induced current density. It follows from il2t that the 
time-periodic part of current density we are interested in does not depend on the permanent part of external current. 
We set /e = again just to simplify the consideration; this is not a limitation of the method employed. 

It is not difficult to find a distribution of surface current density, is{y,t), such that the current is{y,t)[6{x + 
I) — 6{x — I)] shields the superconductor from the external field. Clearly, complete shielding occurs if the magnetic 
vector potential of this current. As, compensates the external magnetic potential inside the superconductor. 



Ae + As = 
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Fig. 3. Standing wave. Tlie external magnetic field is induced by the current Je = ie sin{ky) sm{2nf t)[S{x — a) — S{x + a)]ez. Here 
I = l/fe, L = 10/fc, a = ie = Jc/fc. Shown for t = [0.24, 0.44, 0.56, 1.00, 1.12]// are: the boundaries of plus- and minus 

critical current zones, the magnetic field lines, and the cylinder cross section. Regular finite element mesh 50 x 500, time step 0.04//. 

for —l<x<l. Using expressions il4t and il5l . we find that the shielding would be achieved if we set 

is = -e-'^^^-^he sm{k{y - vt)}, (16) 
is = -e-^("-')ie sin(/cy) sin(27r/t) (17) 

for the running and standing waves, (fTOl i and (fTTT i. respectively. 

We now assume that the depth A to which fluctuations penetrate into the superconductor is much smaller than the 
fluctuation wavelength, kA <^ 1, interpret the shielding surface currents as the integrals of bulk current densities 
across a narrow penetration zone, and find the asymptotic distribution of bulk current density inside this zone 
analytically. First, let us note that for each y the surface current is{y, t) reaches its extremal values when the whole 
penetration zone (see Figs. |2l and |3ll is occupied by the critical current density of the same sign. Therefore, the 
penetration depth can be calculated as A = max(is)/ Jc, which gives A = e-^^^-^^e/Jc for the running wave and 
A(y) = e~'^'("~')ie| sin(A;y)|/ Jc for the standing wave. We see that for both wave types the fluctuations may be 
regarded small if 

z^ = A:Ao<l, (18) 

where Aq = iee-'^^^-'V-^c- 

If, at time t, the current is{y,t) is neither maximal nor minimal, the penetration zone —l<x< — / + A 
(where A = Aq for the running wave and A = Ao|sin(A;y)| for the standing wave) contains regions of plus- 
and minus critical current densities. The near-surface region —l<x<—l + a{y,t) appears at the time when 
dtis{y,t) = and, as it propagates inside, the current density there is JcSign{(9t«s(?/i The rest of the penetration 
zone, — / + (j{y,t) < x < — / + A, is occupied by the critical current density of the opposite sign. Comparing ig 
with the integral of current density across the penetration zone, we find the moving boundary a{y,t): 

a{y,t) = i (^A + ^iMsign{5ii,(y,t)}) . (19) 

Near the slab surface x = +1 the current density distribution is antisymmetric. 

Simple physical arguments were used above to obtain the asymptotic solutions for weak penetration: postulating 
the solution structure, we spread the shielding surface current into a bulk current in the near-surface zone. In a 
similar way, weak penetration of alternating uniform field into a perpendicular circular cylinder has been studied in 
[16], [17]. Although the chosen surface current would have shield the superconductor from the external magnetic 



field completely, spreading this current into the bulk makes shielding imperfect. As will be shown below, the 
remaining field is of the order 0{u). We will now extend our arguments and present the obtained asymptotic 
solutions as zero order terms of consistent asymptotic expansions. 



A. Running wave 

It is not difficult to see that the asymptotic distribution of current density inside the penetration zone —l<x< 
—I + Aq, obtained for u = fcAo ^ 1, can be presented as 

J=-JcS{k{y-vt)-^o{c;/Ao)), <^ = x + /G [0,Ao], (20) 

where S{z) = sign(sin(z)) is a 27r-periodic step-function and "^o{u) = arcsin(l — 2n). 
Let us now look for the current density 

J = -JcSik{y-vt)-^{<;/A)), ? = x + /g[0,A] (21) 

where 

A = Ao(l + riz/ + r2Z/2 + ...), (22) 

^{u) = + z^^'i(u) + iy'^'^2(u) + ... (23) 

are such that the current (l2n i. jointly with the opposite one near another superconductor surface, shield the external 
magnetic field. 

We partly solve this problem in Appendix I by showing first that the vector potential produced by the zero-order 
approximation d20l i to current density compensates the potential of external current inside the superconductor (but 
outside the penetration zone) up to the 0{v) terms; hence the external magnetic field inside the superconductor is 
shielded up to the same order in v. We find then the first order corrections, ri and ^'i, ensuring shielding of the 
magnetic field up to the second order, 0(z^^). 

It turns out (see Appendix I) that a nonsingular function can be found only if ri = — i, so the penetration 
depth becomes 

A = Ao(l-^zv + 0(z.2)). 
It is further shown that the external magnetic field is shielded up to the second order in u if 

M 

^l{c^/A) = «m^2,n(l - 2?/A), (24) 

m=0 

where Unir) = sin((n + 1) arccos(r))/ sin(arccos(T)) are the Chebyshev polynomials of the second kind, 

"™ " {iim + 1)2 - 1 + 4m2 - 1 ) ' ^^^^ 

and M satisfies the conditions 

z^ln(M+ 1) < TT 

if M > 1 (or < vr if M = 0, 1) and 

g-(2M+3)j^ 



< 0(u). 

The first of these conditions ensures that <^ 1, the second provides for the 0{u^) shielding of magnetic field. 

The conditions are easy to satisfy and, usually, only a few terms of the series are needed. Thus, in the examples 
in Fig. H we used Af = 1 for = 0.2 and M = for = 0.4. 



Fig. 4. Penetration of tlie running wave fluctuations. Domains of plus- and minus critical current densities (numerical solution) and tiieir 
asymptotic boundaries: "- -" - zero order approximation, " — " - first order approximation. Left: ie = 0.25Jc/k and u — 0.2; right: 
ie ~ 0.5Jc/k and u = 0.4. In both cases, I — 0.5/fc, a — 0.7/k, L — 2Tv/k. Different scales in x and y. 



B. Standing wave 

We find first a correction to the asymptotic penetration depth A(y) = Aq] sin(A;y)|. Let us choose a moment 
when the external field is the strongest, e.g., t = T/4, so that for each y the induced current density is 

J = -JcS{ky), (26) 

in the whole penetration zone < q = x + I < A{y) (and is opposite in the zone near another slab surface). We 
will now assume that 

A(y) = Ao(| sm{ky)\ + + u^'^2{ky) + ...) (27) 

and find a correction that ensures better shielding of the external field for this moment of time. 

It can be shown (see Appendix II) that the magnetic field is shielded up to the second order in u if 

^i(s) = --sign(sin(s))[2sin(s) +sin(2s)ln(|tan(s/2)|)]. (28) 

TT 

Equations (l26t . (l27l . and d28b give the asymptotic distribution of current for t = T/4. It is now easy to obtain the 
solution for any time moment. Let, for example, T/4 < t < 3T/4. Then the closest to the surface x = —I part of the 
penetration zone is occupied by the current density +JcS{ky). We can present this as a superposition of the current 
density —JcS{ky) in the whole peneti'ation zone and the current +2JcS{ky) in the part of this zone neai^ the surface. 
Up to the second order terms, the current —JcS{ky) shields the external current sm.{ky)[6{x+a)—5{x — a)\, so the 
opposite current density +2JcS{ky) has to shield the external current ie(sin(27r/t) — 1) sin{ky)[5{x + a) — 5{x — a)]. 
This means that the boundary between the two zones must be 

A(y,t) = ^(l-sin(27r/t)). (29) 



In Fig. |5l we compare the numerical and asymptotic solutions for two values of parameter v. As in the running 
wave case, these solutions are close. 




V. AC LOSSES 

Let the external magnetic field be periodic in time, e.g., induced by a periodic current density Je(?^, t), and T its 
period. Then there establishes a periodic induced current density in the superconductor, J(r, t). Suppose this latter 
function was found and it is needed to calculate the energy losses. 



P = ^ [ / J • E dndt = ^ f / pJ^ dndt. 
T Jo Jn ^ Jo Jn 



This would be an easy task for a usual conductor with the known resistivity. However, for type-II superconductors 
p{r, t) is an effective resistivity caused by the movement of magnetic vortices and is not known a priory. Mathe- 
matically, in the Bean model, this is a dual variable excluded by transition to the variational formulation; accurate 
calculation of p (or E) in the general case is difficult. To avoid this complication, we employ a method based on 
the magnetic potential representation of electric field Q and similar to those suggested in [17], [18] but applicable 
to nonuniform external fields. 

Since J satisfies the zero divergence conditions for any gauge holds 

rj^ rj-\ rj~t rj~\ 

P = ^ [ [ 3 ■'Edndt = -^ [ [ J ■dtAdndt = [ [ J-G*dtJdndt-^ [ [ J-dtAedndt. 
T Jo Jn T Jq Jq^ T Jq Jq T Jq Jq 

It is easy to see that 

rT r rT ^ 



n 



dt 2 



J • G * a J dndt = —i- 3-G*Jdn}dt = 



due to periodicity of J. Hence we obtain 

1 '■^ 
/o Jn 

and, since the time period of the product J • dtA.(. is T /2, we can also write 



P = -- I I J-dtAecindt 



2 rto+T/2 

P = -- I 3-dtAedndt (30) 



T 

for arbitrary time moment to. 

A. AC losses for running wave 

Clearly, only half of the slab may be considered due to symmetry and we can now use i30i with = (— ^, 0) 
and T = 2-K/kv to find the asymptotic AC losses for small fluctuations. Obviously, the value of integral (l30t must 
be the same for all y and equal to the rate of AC losses per unit of slab surface. We neglect the second order terms 
in i22l . d23b and assume 

* = ^-o + z.^i, A = Ao(l-^). (31) 
To simplify computations, let us take y = 0, so that (I2H gives 



J\y=o = JcS{kvt + ^(?/A)), ? = X + / G [0, A], 

and choose to = —'^{0)/kv. For small u the function A) is close to ^'o(?/A) and monotonically decreases 
for ? G [0, A]. Therefore, kvto + ^(?/A) changes monotonically from zero at ? = to kvto + ^'(1) = —it at 
<j = A. Indeed, since ^'i(O) = ^i(l), we have 

^(1) = ^o(l) + i^^i(l) = -vr + ^o(O) + i^^'i(O) = -Tr + ^{0). (32) 

Hence, for y = the whole penetration zone at t = to is occupied by the current density —Jc- For to < t < to + T/2 
there is a + Jc-current-density zone propagating inside and sweeping out the — Jc zone at t = to + T/2. The moving 
boundary between the two zones is determined by the condition kvt + ^(?/A) = 0. Taking this into account we 
rewrite d30b as 



2 7 /■'^ ( />->I'(f/A)/fcD i—^{l)/kv 

P = -^ {- dtAedt+ dtAedt 

/o I J-^{0)/kv J-^{i/A)/kv 



T 

2fl0leJC^_^-ka P g-j^ ^(^^^) _ ^^^-^ _ g-j^ ^^^-^^ ^^^^^ ^ ^'^^ 



X — — ; + c; 



kT 



since Ae{x,y,t) is given by (fT4l) . Note that sin^'(l) + sin^'(O) = because of d32l : also ieC ^'^ = JcAqc '^^ 
Substituting A and from d3Tl we obtain 



4^xoJ|Ao 



where, up to the higher order terms. 



Po = -ue ^1 - J sin{arcsin(l - 2u) + u'^i{u)} sinh(-A;/ + uu)du 
^cosh(-/cZ) +sinh(-/c/) j l-i (— ^) a/i - r'^dr 



'^-ki 
2 



(33) 



Here r = 1 — 2u and \l'i is given by 9241 as a finite series of the Chebyshev polynomials. We see that the leading 
AC loss term depends on the first order correction to solution i20l . Since 

' U2mir)Vl^dr = \ ^ "' = JJ' 
I 1^ m > 



only the first term of the series (I24t gives input into AC losses. We find 

,,2 



Po = —e 



-kl 



— cosh(— /c/) ^sinh(— A;/) 

3 2 



and, finally, the asymptotic losses for weak penetration. 



P 



(34) 



where / = l/T is the frequency of fluctuations. Note that the expression obtained would coincide with the well 
known formula for hysteresis losses caused by fluctuations of spatially uniform external magnetic field (incomplete 
penetration), P = "^f ijlqHq/ Jc, if we rewrite this latter formula not for the fluctuation amplitude Hq but using the 
penetration depth Aq, equal to Hq/Jc for the uniform field. 



B. AC losses for standing wave 

To estimate the average asymptotic AC losses per unit of slab surface we modify slightly the formula (I30t to 
calculate the density of these losses, 

P{y)=-- / dtA^Jdat, 

J Jt/4 Jo 

where T = 1//, then average over half the fluctuation wavelength, 

A/2 



{P) = -j^ P{y)dy, 



where A = 2n/k. Assuming T/4 < t < 3T/A and < y < n/k we get J = +Jc for < ? < A{y,t) and 
J = -Jc for A{y,t) < A(y), where A{y) = Ao(| sin(A;y)| + V^i{ky) +0{v'^)) and A(y,t) is given by d29l . 
Substituting 1^ = ^, x = —I + <; into the vector potential (ITsl and integrating, we obtain 



2 J r'^^i'^ 
P{y) = -4^ / 



-A(s/,t) 



where 



^o(y) 



-kl 



sin(/cy) 



cosh(— /cZ) + cosh ( —kl + v 



A(t/,t) 



.A(y) 



dtAedii \ dt = i^fiofJ^AoPoiy), 



An 



sinh (—kl + i^-^^ 



s'mh(—kl) 



2Ao 



(35) 



For < y < ir/k we have A(y)/Ao = sin(A;y) + zv^i(A;2/) + 0{v'^) and, up to the higher order terms. 



^'o(y) 



-kl 



— sin (/cy) cosh(/i;/) + ( — ^i(/cy)sin {ky) cosh.{kl) 



1 

12 



sin^(A;y) sinh(A;/) 



Except for the points very close to y = or vr/A;, where the term containing ^i(A;y) dominates but the losses are 
negligible, 

-2kl 

— • (36) 



P{y) « ^Aio/Jc'Ao^'sin3(A;y)e-^'cosh(A:/) = l^^f Jl^l{y)s\n\ky)'^^ 



This is close to the density of losses in an alternating uniform field ii kl <^ 1 and is twice smaller if A;Z » 1, 
provided the penetration depth is the same. Finally, we find the asymptotic average AC losses per unit of slab 
surface for the standing wave fluctuations: 



(37) 



We see that in this case the leading asymptotic AC loss term is determined by the zero order approximation to 
current density distribution and does not depend on the first order correction as in the running wave case. 



VI. Discussion 

Magnetic field fluctuations are inevitable in most practical applications of superconductors. In this work we used 
the Bean critical state model to study the effect of spatially nonuniform fluctuations of the external magnetic field. 
Although solutions have been obtained only for two model situations, they help to understand qualitatively the effect 
of nonuniform fluctuations in general. The asymptotic estimates, derived for small spatially nonuniform fluctuations 
using a general gauge-invariant formula for AC losses, are the main result of our work. 

For an alternating uniform magnetic field the AC losses in a superconductor are usually expressed via the ampli- 
tude of field variations. This is inconvenient if the external field is nonuniform. Typically, as for the configurations 
considered in our work, both the magnitude and direction of the external field depend on position, time, problem 
geometry, and the fluctuation wavelength. It seems difficult to relate the losses to any specific characteristic of this 
field. The lack of a universal direct relation between the external field at the surface of a superconductor and AC 
losses becomes even more apparent if we consider, for example, a hollow superconductive cylinder with a long coil 
placed into its hole [19]. An alternating coil current induces the shielding current in the superconductor because the 
magnetic flux changes. Nevertheless, had the superconductor been removed, the same coil current would produce 
no magnetic field outside the coil at all. 

Because of this reason the formulas for AC losses in this work are presented in terms of the depth to which 
fluctuations penetrate into the superconductor. To determine this depth and the induced current density asymptotically 
for small fluctuations, we found first the shielding surface current. Spreading this current into the bulk and taking 
into account the current density constraint, we were able to obtain the zero order approximation to current density 
which was then further improved. 

Such an approach is not limited to slab configuration and the two types of fluctuations considered above. Thus, 
using the method derived recently by Bhagwat and Karmakar [20], one can find the surface current shielding a 
cylinder of an arbitrary given cross section placed into a uniform external magnetic field. This makes possible 
to extend, following the scheme used in our work, the asymptotic solution for cylinders in alternating uniform 
transverse field [16], [17] to cylinders with non-circular cross sections (in this case the penetration is weak if its 
depth is much smaller than the characteristic cross section size). 

The results obtained enable one to estimate AC losses in superconductors of magnetic bearings and levitation 
systems, where the typical configuration is similar to that of a running wave fluctuations of external field near the 
surface of a thick slab {kl » 1, one-sided action of a nonuniform external field). Suppose we can remove the 
superconductor and measure the tangential component of field fluctuations at the position of slab surface. Let us 
approximate this component by a running wave with an amplitude Ht. By the method of images, surface current 
shielding these fluctuations has the amplitude 2Ht. This can be used to estimate the penetration depth, Aq ~ 2Ht/ Jc- 
The losses can now be approximated using the formula (I34t derived for the running wave fluctuations. 

Numerical simulations based on a variational reformulation of the critical-state model helped us to envision 
solution structures and to control accuracy of asymptotic solutions. The asymptotic solutions, obtained at first by 
means of simple physical arguments, were presented as zero-order terms of a consistent asymptotic expansion. 
Finding the first order correction allowed us to improve these solutions. It has been shown that the correction 
ensures shielding of external magnetic field up to the second order, H/Hf, ~ 0(i^^), and (see Figs. |4l|5ll provides 
for a satisfactory approximation for small parameter values up to ~ 0.4. For both the running and standing wave 
fluctuations, the maximal penetration depth Amax is smaller than the value Aq given by zero order approximation: 
we found /S.max ~ Ao(l — f) in the first case, A^ax ~ ^o(l — v) second one. 

Expressed via the penetration depth, AC losses for running wave fluctuations d34l l are given by exactly the same 
formula as AC losses for uniform field fluctuations. We cannot expect such a coincidence also for standing wave 
fluctuations, since the penetration zone depth is not constant in that case. However, the local losses (I36t are close 
if /c/ <C 1, which simply means that, locally, the long- wave limit is similar to the uniform fluctuations case. If 
the wavelength is shorter, the losses can be at most twice smaller than in the uniform case for the same local 
penetration depth. Using the average penetration depth, (A) = 2Ao/vr, we can present the average losses (l37l as 
(P) = |C/io/Jc (A)3, where C = 7r2(l + e-'^^^)/12. Since 0.82 < C < 1.64, the average losses are also close to 
AC losses in the uniform field case. 

For standing wave fluctuations, the leading AC loss term is determined by zero order approximation to current 
density distribution, i.e., an approximation that is comparatively easy to find. This is similar to the case of uniform 



field fluctuations considered earlier in [17]. It is, therefore, surprising that the situation is different for the running- 
wave type of fluctuations: here knowledge of the first order correction to such a solution is necessary. Neglecting 
this correction would cause an error of the same order as the AC loss value itself. 

To explain this discrepancy, let us note that the expressions i33l and d35l) for AC losses for running and 
standing waves, respectively, have different structures. In d35l) . the small parameter appears only in a combination 
z/A(y)/Ao = i^^'(y), where = + z^^^i + ... We can write Po{y) = J^{y, ^'^{y)) and check that the function 
T{y,u) satisfies J^{y,0) = J^u{y,0) = 0. Hence Po{y) = ^J^uu{y,0)^o{y) + 0{v^); the main term does not 
depend on ^fi. 

In the running wave case, the dependance on small parameter is different because the penetration depth and 
the shape of the free boundary are separated. The equation (l33l can be written as Pq = vT{y^ ^'), where is a 
functional and ^' = ^'o + ^^^i + ••• as in the previous case. Expanding we get Pq = ^'o) + v^\Tv(^^ ^'o) + 

(.F.i>(0, ^'o), ^i)] + 0{y^). Here Ti^ is the Frechet derivative, the first term turns out to be zero, and so the leading 
term depends on ^i. 

In both cases, however, the leading term of AC losses is proportional to Aq, which corresponds to P 
dependance known for the uniform field fluctuations. The next approximation would lead to a deviation from cubic 
law. In the frame of the Bean model, such deviation can be due to the shape of a superconductor being different from 
that of a slab [6], [19], heating caused by AC losses [6], or, as in the present case, because of spatial non-uniformity 
of the external field. 

Appendix I 

The vector potential of the current dSTT i and the opposite one near x = +/ can be written as 

^. = . J = £ £ - .t) - *(./A)) in (fe±i^|±|5^) * 

Changing the variables, u = ?/A, s = k{y' — vt) — ^(n), and using the Taylor expansion we obtain: 



A, = ^Ao(l + nv + ...) r S{s) I In 



\ [k{x -l) + v^uf + [k{y -vt)-s- *(n)]2 




2k'n: 
Here 



[Ao + Aiu + 0{iy^)] . (38) 



ho = |e-^('^-^) (39) 
is introduced as a characteristic magnitude of external field fluctuations at the surface of the superconductor and 

Aq = D{x,y,t,u)du, Ai = uFi{x,y,t,u)du + / '$i{u)F2{x,y,t,u)du + viAq, 
Jo Jo Jo 



where 



77 _ _o / q( \ I + - , , 

^" / ^'^\iT-sy + [kix + l)]^^{T-s)^ + [kix-l)]\l ' 



oo 



r = k{y-vt) -^o{u). 
To calculate the integrals we present the periodic function S as the Fourier series. 



oo ^ 

S{s) =^bnsm{ns), bn = I 

n=l ^ 



— n = 2m + 1, 
n = 2m 



and obtain 



oo 



D = 47r ^ ^g-nfci ginh(nA:x) sin(nr), 

n=l 

oo oo 

Fi = 47r ^ hnC""'"' sinh(n/cx) sin(nr), F2 = -47r ^ ftne""''' sinh(nA;x) cos(nr). 

n=l n=l 

We further obtain 

00 ^ 

= 47r — e""^' sinh(nA;j;)[a„ sin(nA:(y — ut)) — /3„cos(n/c(y — ft))], 
n 

n=l 

where 6n = for all even n and, for odd n, 

/■I C 71-/4 n = 1 

an = cos{n"^Q{u))du = < ' Pn = sin(n^'o(u))dn = 0. 

Jo [ U n > i 



Hence, 
Similarly, 



Ao = 47re sinh(A:a;) sin(/c(y - vt)). (40) 



Ai =47r^^ 6^6 "'^^ sinh(nA;x) [7„ sm{nk{y — vt)) — 6n cos{nk{y — vt))] 



n=l 



(41) 



+ ri47re '''' sinh{kx) sm{k{y — vt)) , 
where all even series terms are zero and, for n odd. 



7n = / ucos{n"^Q{u))du — / ^i(n) sin(n'I'o(n))du, (42) 
Jo Jo 

6n= usm{n'^o{u))du + / ^'i(n) cos(n^'o(ii))(in. (43) 
Jo Jo 

We can now substitute the expressions (l40l and d^TT i into (l38l . use (I14t with 1^ = 0, and calculate the total magnetic 
field inside the superconductor but outside the penetration zone. We see that A = A^ + Ai = ^^^(j^^i + 0{i''^)). 
Therefore, up to the second order terms, 

Hrc =2z^/io{rie"^' sinh(A;x) cos{k{y - vt))+ 



nbne "'^^ sinh(n/cx) [7„ cos{nk{y — vt)) + 6n sm{nk{y — vt))]}, 

n=l 

Hy = - 2uhQ{rie~^'- cos\i{kx) sin{k{y - vt)) + 
00 

nfe^e""^^ cosh(nA;x)[7„ sin(nA;(y — vt)) — 5nCos{nk{y — vt))]}. 

n=l 

This proves that the current density i20l shields the external field up to the first order in v and is a zero order 
approximation. To nullify the first order terms of magnetic field we will now try to satisfy the conditions 

— ri/61 n = 1, , ^ ^ .... 

n>i; ^"'^ ^- = 0' (44) 

for odd values of n. Since ^'o(^) = arcsin(l — 2n), we denote r = l — 2u and express sin(n^'o(^)) and cos(?7,'I'o(^^)) 
in (I42t . (l43l via the Chebyshev polynomials of the first and second kind, 

Tn{T) = cos(n arccos r), Un^r) = sin((n + 1) arccosr)/ sin(arccos r), 

correspondingly, for r G [—1, 1]. It can be shown that, for n odd, 

sin(narcsinr) = (— l)~T„(r), cos(narcsinr) = {—l)~Un-i{T)\/l — r^. 



Defining ipiiT) ^'1(^-2^) = ^i{u) and calculating the integrals of known functions, we obtain 



7n 



i(-l)^ (dn + J\ Mr)Un-i{r)VT^dT'^ , 



2 

2 



where n is odd and 



dn = l 



+ 



(n + 1)2-1 (n- 1)2-1 



We can now rewrite the conditions (l44l as 



%i(T)r„(r)cir = | 4(^+^1) ^^J' (45) 



1 

1 



Mr)Un-i{rWl-T^dT = -dn (46) 



1 

and use them to determine iIji{t) and ri. Let us present V'l as the sum of its even and odd parts, ipi and ifj". Since 
functions are odd for odd n, even for even n, and orthogonal on [—1, 1] with the weight l/^fl^^T^, condition 
(l45l means that ^°(r)\/l - r2 = cri(T) where Ti(t) = r and 

Thus ^/'j' = cr/\/l — r2 and is singular at r = ±1 if c 7^ 0. Since the expansion (l38l is valid only for |z^'I'i(u)| <C 1, 
we need a nonsingular function and must set 

ri = 

2 

to make c = 0. This determines the first order correction to the penetration depth: 

A = Ao(l-^zv + 0(z.2)). 

The function ipi becomes even and we can expand it into a series of the Chebyshev polynomials of the second 
kind containing only the polynomials of even orders: 

00 

= X] amU2m{T). (47) 
m=0 

The polynomials C/„ are orthogonal on [—1, 1] with the weig ht ^/l — t2, so the coefficients of this series are easily 
found from the condition (I46t : 

= -^^2^+1 / ULir)V^'dr = -i (4(^/^)2 _i + ■ 

We have now satisfied the conditions (l44li but there appears a contradiction: although the series (l47t converges for 
r G (—1,1), t/ji becomes infinite at r = ±1. Indeed, at these points U2m = 2m + 1 and Um ~ — l/27rm2 for big 
m. Since the conditions (I44t were derived under assumption |i/^i(n)| <C 1, the function determined by this series 
cannot be accepted as the first order correction to "^q. 

The computations, however, were not in vain and this singularity can be eliminated, if we take only a finite 
number of terms and define 

M 

V'i(^) = Yl «™f^2m(r). (48) 

m=0 



It can be shown that in this case < ln(M + 1)/tt for all M > 1 and |^i(r)| < 16/157r for M = or 1. 

Thus, the assumption remains valid provided that 

i/ln(M + l)<7r ifM>l, 
1/ « vr if M = 0, 1. 

It is easy to see that conditions (l44l i for 7„ are still satisfied for all n, whereas the conditions for 5„ hold only for 

n < 2M + 1. For n > 2M + 1 we have 

The non-compensated magnetic field inside the superconductor has, up to the second order in v, the components 

oo 

Hx = 2vhQ nfen^nC""^' sinh(n/cx) sin(n/c(y — vt)), 

n=2Af+3 
oo 

Hy = 2vhQ nbn5ne~^^^ cosh(n/cx) cos{nk{y — vt)). 

n=2M+3 

For |x| < / — Aq we have e""'^' cosh(nA;x) < e^"'^, e~^^^\ smh.{nkx)\ < e~^'^; also n6„ = 4/7r for odd n. Hence, 

m=A/+l 



vr (2m + 1)2 % Jj^i (2m + 1)2 %(2M + 1) 



which proves that the field can be made 0{v'^) if 

(2M+3)!. 

This means that the first order correction to '^q{u) may be chosen as '^i{u) = ^^1(1 ~ 2n), where ijji{T) is given 
by (l48l and Af satisfies the conditions d49t and dSOt . 

Appendix II 

The vector potential of the current (l26t can be written as 



Changing the variables, u = ?/A(y') and s = we obtain 

A^nJr r f /[A;(x + /)-i/^^u]2 + (s-A:y)2\ 1 

Ai = ^ S{s)A{s/k) I Ini — -r^^ ^ — ]du}ds 

4fc7ry_^ ^' '\Jo \[k{x-l) + u^u]^ + {s-ky)y J 

with S{s)A{s/k) = Ao(sin(s) + z^^'i(s)S'(s) + ...). Integrating in u and expanding we get 

where the characteristic amplitude of fluctuations /iq is defined by d39b . 

f[k{x + l)]^ + {s-ky)^\ . , _kl ■ n ^ ■ un ^ 

Aq = In -^-^ ;(^75 ) sm(s) ds = 47re sm(A;y) smh(kx), 

J-00 \[k{x - l)]^ + {s - ky)y 

A, = I Ms)s{s) In j 



00 



loo I [k{x + + is- ky)^ + [Mx - /)]2 + (s - %)2 / 



(51) 



Obviously, ^'i should be a vr -periodic and even (because of symmetry) function. To calculate the integrals in dSTl 
we present ^'i(s)S'(s) and | sin(s)| sin(s) as the Fourier series, 



^'i(s)S'(s) = an sin( 



ns 



sin(s)| sin(s) = 6„ sin( 



ns 



n=l 



n=l 



where bn = ~ 7rn(ra^-4) ^'^^ ^'^'^ n, bn = for even n, and the coefficients an are unknown. Integrating, we 



find 



Ai = 27r^ i-an + bnj 

n=l ^ ^ 



-knl 



s'm{nky) sinh(nA;a;) 



Up to 0(z^^), the total magnetic potential can be written as 



zkir k 



n=l ^ ^ 



e '^"^ sin(n/c?/) sinh(nA;x). 



To make zero the 0{v) terms of magnetic field inside the superconductor, it is sufficient to satisfy the conditions 



-n- 



7r(n2-4) 





n = 2m + 1, 
n = 2m. 



This gives 



^'i(s)5(s) = ^ a2m+i sin({2m + l}s) 



m=0 



and so, up to the second order terms, the penetration depth can be presented as 

= Ao{| sm{ky)\ + z^y(A:y)sign[sin(A;y)]}, 

where the series 

4 ^ sin({2m+l}s) 
^ ' TT ^ (2m+l)2-4 

m=0 ^ ^ 

can be summed up analytically. Let us denote z = e**, 

4 ^ cos({2m + l}s) 
^ ' TT ^ (2m + 1)2 -4 ' 

m=0 ^ ^ 

and consider the complex function 

' — ' „ , , oo 



w =u{s) + iv{s) = -y - 

vr ^ (2m + 1)2-4 vr 

m=0 m=0 



1 



1 



1 

vr 

1 

vr 



--Z + { z^ 

z 



1] V ^ 

z2 / ^ 2m - 1 

^ m=l 



1 

vr 



--z+[z^ 

z 



2m - 1 2m + 3 

oo 



^2) 



m=l 







- - Z + ( - 4t 



Z' 

Since 2; — l/z = 2isin(s), — l/2;2 = 2isin(2s), and Re 
W is 



r dz 


1 


"1 


Jo 


vr 


z 



1 \ 1 



1 + z 



z^ 2 \l-z 



- - z + { z^ - ^ ] - In 



ln(| tan(s/2)|), the imaginary part of 



2 1 

- sin(s) sin(2s) ln(| tan(s/2)|). 

vr vr 
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